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Abstract. We study a two-scale reaction-diffusion system with nonlinear re- 
action terms and a nonlinear transmission condition (remotely ressembling 
Henry's law) posed at air-liquid interfaces. We prove the rate of convergence 
of the two-scale Galerkin method proposed in [7] for approximating this sys- 
tem in the case when both the microstructure and macroscropic domain are 
two-dimensional. The main difficulty is created by the presence of a boundary 
nonlinear term entering the transmission condition. Besides using the par- 
ticular two-scale structure of the system, the ingredients of the proof include 
two-scale interpolation-error estimates, an interpolation-trace inequality, and 
improved regularity estimates. 



1. Introduction 

Reaction and transport phenomena in porous media are the governing processes 
in many natural and industrial systems. Not only do these reaction and trans- 
port phenomena occur at different space and time scales, but it is also the porous 
medium itself which is heterogeneous with heterogeneities present at many spatial 
scales. The mathematical challenge in this context is to understand and then con- 
trol the interplay between nonlinear production terms with intrinsic multiple-spatial 
structure and structured transport in porous media. To illustrate this scenario, we 
consider a large domain with randomly distributed heterogeneities where complex 
two-phase-two-component processes are relevant only in a small (local) subdomain. 
This subdomain (which sometimes is refered to as distributed micro structur^ fol- 
lowing the terminology of R. E. Showalter) needs fine resolution as the complex 
processes are governed by small-scale effects. The PDEs used in this particular 
context need to incorporate two distinct spatial scales: a macroscale (for the large 
domain, say SI) and a microscale (for the microstructure, say Y). Usually, x G Q 
and y £ Y denote macro and micro variables. 

1.1. Problem statement. Let S be the time interval ]0, T[ for a given fixed T > 0. 
We consider the following two-spatial-scale PDE system describing the evolution of 
the the vector (U,u,v): 
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(1.1) 6d t U(t,x) -DMJ{t,x) = - / b(U(t,x) -u(t,x,y))d\ 2 y in S x SI, 

(1.2) d t u(t,x,y) - diA y u(t,x,y) = -kr)(u(t,x,y),v(t,x,y)) inSxSlxY, 

(1.3) dtv(t, x, y) — d,2A y v(t, x, y) = — akr){u[t, x, y), v(t, x, y)) in 5 x SI x Y, 

with macroscopic non-homogeneous Dirichlct boundary condition 

(1.4) U(t,x) = U ext (t,x) on SxdSl, 
and microscopic homogeneous Neumann boundary conditions 

(1.5) Vj/tt(t, x, y) ■ n y = on S 1 x J] x Tjv, 

(1.6) V y «(£, x, y) • % = on S x SI x I\ 

The coupling between the micro- and the macro-scale is made by the following 
nonlinear transmission condition on Tr 

(1.7) — S7 y u(t, x, y) ■ n y = — b(U(t, x) — u(t, x, y)) on S x Q x Tr. 
The initial conditions 

(1.8) £7(0,a:) = Ui{x) in O, 

(1.9) u(0, x,y) — Ui(x,y) in fi x Y, 

(1.10) v(0, x, y) = vi(x,y) in CI X Y, 

close the system of mass-balance equations. 

Continuing along the lines of [7] , the central theme of this paper is understanding 
the role of the nonlinear term b(-) in what the a priori and a posteriori error analy- 
ses of (|l.ip - (|1.10p are concerned. Within the frame of this paper, we focus on the a 
priori analysis and consequently prepare a functional framework for the a posteriori 
analysis which is still missing for such situations. Since our problem is new, the 
existing well-established literature on a priori error estimates for linear two-scale 
problems (cf. e.g. [5]) cannot guess the rate of convergence of the Galerkin approx- 
imants to the weak solution to (ll.ll) - (|1.10p . Therefore, a new analysis approach is 
needed. Notice that the main difficulty is created by the presence of a boundary 
nonlinear term entering the transmission condition (jl.7p . Here we prove the rate 
of convergence of the two-scale Galerkin method proposed in [7] for approximating 
this system in the case when both the microstructure and macroscropic domain 
are two-dimensional, see Theorem 13.51 Nevertheless, we expect that the results 
can be extended to the 3D case under stronger assumptions, for instance, on the 
regularity of Tr and data. Besides using the particular two-scale structure of the 
system, the ingredients of the proof include two-scale interpolation-error estimates, 
an interpolation- trace inequality, and improved regularity estimates. 

The paper is structured in the following fashion: 
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1.2. Geometry of the domain. We assume the domains fl and Y to be connected 
in IR 3 with Lipschitz continuous boundaries. We denote by X k the A;-dimensional 
Lebesgue measure (k G {2, 3}), and assume that A 3 (£7) ^ and A 3 (F) ^ 0. Here, 

is the macroscopic domain, while Y denotes the part of a standard pore associated 
with microstructures within f2. To be more precise, Y represents the wet part of 
the pore. The boundary of Y is denoted by T, and consists of two distinct parts 

r = t r u rv 

Here r^nTjv =0, and A^Tr) ^ 0. Note that T N is the part of dY that is isolated 
with respect to transfer of mass (i.e. Tjv is a Neumann boundary), while Tr is 
the gas/liquid interface along which the mass transfer takes place. Throughout the 
paper Xy (k S {1, 2}) denotes the fc-dimensional Lebesgue measure on the boundary 
dY of the microstructure. 

1.3. Physical interpretation of (|1.10|) . U, u, and v are the mass concen- 
trations assigned to the chemical species Ax, Ai, and A$ involved in the reaction 
mechanism 

(1.11) Ax ±=r A 2 + A 3 H 2 + products. 

For instance, the natural carbonation of stone follows the mechanism (jl.lip , where 
Ai := C02(g), A 2 := C02(aq), and A 3 := Ca(OH) 2 (aq), while the product of re- 
action is in this case CaC03(aq). We refer the reader to Q] for details on the 
mathematical analyis of a (macroscopic) reaction-diffusion system with free bound- 
ary describing the evolution of (jl.llj) in concrete. 

Besides overlooking what happens with the produced CaC03(aq), the PDE sys- 
tem also indicates that we completely neglect the water as reaction product in 
(|1.11|) as well as its motion inside the microstructure Y. A correct modeling of the 
role of water is possible. However, such an extension of the model would essentially 
complicate the structure of the PDE system and would bring us away from our 
initial goal. On the other hand, it is important to observe that the sink/source 
term 

(1.12) - f b{U-u)d\ 2 y 

models the contribution in the effective equation (|1.1[) coming from mass transfer 
between air and water regions at microscopic level. Surface integral terms like (|1.12j) 
have been obtained in the context of two-scale models (for the so-called Henry and 
Raoult laws [3] - linear choices of b(-)l) by various authors; see for instance [5] 
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and references cited therein. The parameter k is the reaction constant for the 
competitive reaction between the species A2 and A3, while a is the ratio of the 
molecular weights of these two species. Furthermore, we denote by 9 the porosity 
of the medium. 

2. Technical preliminaries 

2.1. Assumptions on data, parameters, and spatial domains Q,,Y. For the 
transport coefficients, we assume that 

(Al) D > 0,di > 0,d 2 > 0. 
Concerning the micro-macro transfer and the reaction terms, we suppose: 

(A2) The sink/source term b : M — >• M + is globally Lipschitz, and b(z) = if 

z < 0. This implies that it exists a constant c > such that b(z) < cz if 

z > 0; 

(A3) r] : R x R -> M+ is defined by jj(r,s) := R(r)Q(s), where i?,Q are glob- 
ally Lipschitz continuous, with Lipschitz constants cr and cq respectively. 
Furthermore, we assume that R{r) > if r > and i?(r) = if r < 0, and 
similarly, Q(s) > if s > and Q(s) = if s < 0. 
Finally, we have k > 0, and a > 0. 
For the initial and boundary functions, we assume 

(A4) u ext g ii 1 (5, J H" 2 (n))nif 2 (S',i> 2 (n))n£y(»s x n), c/> g f 2 (0) ni^(O), 

17/ - U ext (0, •) G T/o 1 ^), uj, G 7 2 (0, i7 2 (F)) n L^°(0 x F). 
For the approximation with piecewise linear functions (finite elements), we assume: 

(A5) f2 and Y are convex domains in R 2 with sufficiently smooth boundaries; 

(A6) h 2 max{7i, 73} < 1, where and 73 are strictly positive constants en- 

tering the statement of Lemma 13.11 

2.2. Weak formulation. Known results. Our concept of weak solution is given 
in the following. 

Definition 2.1. A triplet of functions (U,u,v) with (U - U ext ) E L 2 (S,H^(n)), 
d t U G L 2 {S x Q), (u,v) G L 2 (S,L 2 (n,H 1 (Y))) 2 , (d t u,d t v) € L 2 (S xtlx Y) 2 , is 
called a weak solution of (|l.l[l - (|1.10[) if for a.e. t £ S the following identities hold 

(2.1) — / 6Up+ [ D\7UVip+ [ [ b(U - u)ipd\ 2 y dx = 

d 
di 

(2.2) +k I r)(u,v)(t> = 

JnxY 

(2.3) ^- [ vip+ [ d 2 VyvV y ip + ak [ j](u,v)ip = 0, 
dt JnxY JnxY JnxY 

for all (</>,</>, V) G ^o( n ) x l2 {^ h1 ( y )Y, and 

[/(0) = Uj in u(0) = u/, u(0) = «/ in [1x7. 

Theorem 2.2. 7i exists a globally-in-time unique positive and essentially bounded 
solution (U,u,v) in the sense od DeHnition \2. 1[ 

Proof. We refer the reader to [7] for the proof of this result. □ 



si Jn JnJr a 

2 



u4>+ diV y uVy<j)- / / b(U - u)(j)d\ y dx 

nxY JnxY JnJr R 
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2.3. Galerkin approximation. Basic (semi-discrete) estimates. Following 
the lines of [319], we introduce the Schauder bases: Let {£i}ieN be a basis of L 2 (fl), 
with £j £ -ffg(f2), forming an orthonormal system (say o.n.s.) with respect to 
L 2 (S!)-norm. Furthermore, let {C.j'fc}i,fceN be a basis of L 2 (VL x Y), with 

(2-4) Cjk(x,y) = £j(%)r)k(y), 

where {r/k}k&i is a basis of L 2 (Y), with r/k 6 H 1 ^), forming an o.n.s. with respect 
to L 2 (Y)-norm. 



5/V pN 

associated to the bases {CjljeN, and {r)k, respectively. For (y>, ip) of the form 



Let us also define the projection operators on finite dimensional subspaces P£ , Py 



j.keti 



we define 

(2.5) = X>&( a 



iV 



(2.6) (P?1>)(x,v) = DEiiK7iWi»(y) 



JV 



(2.7) {p?i>)(*,v) = EE^w*)^)- 

jeN fe=i 



The bases {<7j}jeN> and {?7fe}fcgN are chosen such that the projection operators 

oN pN 



P^ \Py are stable with respect to the L°°-norm and £f 2 -norm; i.e. for a given 



function the L°°-norm and i? 2 -norm of the truncations by the projection operators 
can be estimated by the corresponding norms of the function. 

Remark 2.3. Apparently, this choice of bases is rather restrictive. It is worth noting 
that we can remove the requirement that ,Py are stable with respect to the 
i°°-norm in the case we work with a globally Lipschtz choice for the mass-transfer 
term &(•). We will give detailed explanations on this aspect elesewhere. 

Now, we look for finite-dimensional approximations of order N G N for the 
functions Uq :— U — U ext ,u, and v, of the following form 



(2.8) U N (t,x) = £ of («&(*), 

TV 

(2-9) u N (t,x,y) = E (3?k(m(x)Vk(y), 



i,k=i 

N 



(2-10) v N {t,x,y) = E ^m^)v k (y), 

J,k=l 
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where the coefficients , [3^ k , 7^, j, k = 1, . . . , N are determined by the following 
relations: 

(2.11) / 9d t U^ (t)ipdx + j D\7U^{t)\7ipdx = 



yj b ((Uq + U ext - u N )(t)) d\ 2 y + 0d t U ext (t) + DAU ext (t)^J tpdx 



in \Jr f 

(2.12) / d t u N {t)(j)dxdy + / d^yu" (t)\7 y <j)dxdy 



nxY JnxY 

b ((C/ W + U ext - u N ){t)) <\>d\\dx - k f t) (u N (t),v N (t))(f)dyda 

lnJT R JnxY 

(2.13) / d t v N \t)ip dydx + [ d 2 \7 y v N '(t)V„V dydx = 
JnxY JnxY 

- ak f n(u N (t),v N (t))^dydx 
JnxY 



for all if G spanj^j : 


3 G {1, 


and 




(2.14) 


<(0) 


(2.15) 


/$(0) 


(2.16) 


TjX(O) 



(tfj - (7 e **(0))^x, 
ui(j k dxdy, 
viQkdxdy. 



n jy 



>n jy 

Theorem 2.4. Assume that the projection operators ,P y , defined in 
\2. 7| ), are stable with respect to the L°°-norm and H 2 -norm, and that (Al)-(A^) 
are satisfied. Then the following statements hold: 

(i) The finite- dimensional approximations U^(t), u N (t), and v N (t) are posi- 
tive and uniformly bounded. More precisely, we have for a.e. (x,y) S fixY , 
all teS, and all N E N 

(2.17) < U^(t,x) < mi, <u N (t,x,y) < m 2 , < v N (t,x,y) < m 3 , 

mi := 2||[/ e;rt || L oo( Sx0) + ||C//|| L oc(n), 
m 2 := max{||u/|| i oc( f2x y ) ,mi}, 
m3 := ||«/|U~(nxK)- 

(ii) There exists a constant c > 0, independent of N, such that 

( 2 - 18 ) l|C / o V ||i°°(s,i/ 1 (n)) + \\dtUo\\ L 2( SiL 2 {n)) < c, 

(2.19) ||u A '||l^(S,L 2 (0;H 1 (^))) + \ \9tU N \ \L2(S,L 2 (n-,L 2 (Y))) < c, 

(2.20) ||« Ar ||L°°(S,L 2 (n ; ff 1 (y))) + ||&« IU»(,^£ a («;£ 2 (y))) < c ; 

(hi) Then there exists a constant c > 0, independent of N, such that the follow- 
ing estimates hold 

(2-21) W^xU N \\ L ^(^s,L 2 (nxY) + ||Va:V ||x,»o(S ,L 2 (f2xF) < c 

(2.22) ||VyV x U Ar || L 2( S i 2( 0)< y) + ||VyV X W Ar || L 2( S L 2 (nx y) < C. 
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Proof. This statement combines the information stated in Theorem 6.1 and Theo- 
rem 6.2 from [7j. We refer the reader to the cited paper for the proof details. □ 

With these estimates in hand, we have enough compactness to establish the 
convergence of the Galerkin approximates to the weak solution of our problem. 

Theorem 2.5. There exists a subsequence, again denoted by (Uq , u N ,v N ), and a 
limit (U ,u,v) G L 2 (S;H 1 (n))x[L 2 (S;L 2 (n-,H 1 (Y)))] 2 ! with(dtU^,d t u N ,d t v N ) G 
L 2 (S xfl)x [L 2 (S xflx Y)] 2 , such that 

(U^,u N ,v N ) -> (U ,u,v) weakly in L 2 (S '; H 1 (ft)) x [L 2 (S; L 2 {Vl- H l {Y)))] 2 
(dtl/Q , d t u N , d t v N ) —> (dtUo,dtu,dtv) weakly in L 2 
(Uq ,u , v N ) (Uo,u,v) strongly in L 2 
u N \r R — > u\r R strongly in L 2 (S x fl, L 2 (r#)) 

Proof. See the proof of Theorem 6.3 in [7J. □ 

In the next section, we address the question we wish to answer: 

How fast do the subsequences mentioned in Theorem \2.4\ converge to their unique 
limit indicated in Theorem \2.,5\ ? 

3. Estimating the rate of convergence: The case Y c ft c M. 2 

Adapting some of the working ideas mentioned in [TU1 H] to this two-spatial-scale 
scenario, we obtain an a priori estimate for the convergence rate of the Galerkin 
scheme constructed in section [2~3l 

3.1. Approximation of smooth two-scale functions. As preparation for the 
definition of the finite element solution to our problem, we briefly introduce some 
concepts concerning the approximation of smooth functions in f2, Y C 1R 2 (taking 
into account assumption (A5)); see, for instance, [2] or [TU] for more details. 

For simplicity, we let h denote the maximum length of the sides of the triangu- 
lations Th of both 17 and Y . h decreases as triangulations are made finer. Let's 
assume that we can construct quasiuniform triangulations ([ID], p. 2) and that the 
angles of these triangulations are bounded from below by uniformly in h positive 
constants. 

Define Vh ■= span{£j : j E {1, . . . , N}}, and Bh '•= spanf^ : k G {1, . . . , N}} 
where and r]k are defined as in section [2~3l We also introduce Wh '■= span{Cjfe : 
j, k G {1, . . . , N}}, where ( jk are given by (f!Q|) . Note that W h ■= V h x B h . 

A given smooth function ip in Q vanishing on dfl may be approximated by 
the interpolant l^Lp in the space of piecewise continuous linear functions vanishing 
outside (J Th ■ Standard interpolation error arguments ensure that for any ip G 



\\V(I h (p - ^)||i2(n) < ch\\(p\\ L 2 {n} . 
We define the macro and micro-macro Riesz projection operators (i.e. IZff and 
72.™) in the following manner: 



H 2 (n)nH^ 



(ft), we get 



Whip - ip\\ L 2 (n) < c/i 2 ||v3|| L 2 (0) 



(3.1) 
(3.2) 




K% :L 2 {n;H\Y))^W h , 
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where R¥ is the standard single-scale Riesz projection, while is the tensor 
product of the projection operators 

(3.3) P m : L 2 (n)^V h 

(3.4) P n : H^iY) -)■ B h . 



Note that this construction of the micro-macro Riesz projection is quite similar 
to the one proposed in [6] (cf. especially the proof of Lemma 3.1 loc. cit.). The only 
difference is that we do not require any periodic distribution of the microstructure 
Y. Consequently, if one assumes a periodic covering of f2 by replicates of Y sets, 
then one recovers the situation dealt with in [BJ. 

Lemma 3.1. (Interpolation-error estimates) Let 1Z™ andlZff be the micro and, re- 
spectively, macro Riesz' s projection operators. Then there exist the strictly positive 
constants 7^ (i G {1,2,3,4}J, which are independent of h, such that the Lagrange 
intepolants 1Z™4> and U-ff ip satisfy the inequalities: 



(3-5) \\<p-Kh<p\\L 2 (n) < 7ifr 2 |Mlff 2 (n), 

(3-6) \\tp -Tl^(p\\ H i (n) < j2h\\(p\\ H 2(n), 

(3.7) \\(f - 1?h<l>\\L*(n;L*{Y)) < 73^ 2 (\\4>\\L 2 (n:H 2 (Y))nL 2 (Y;H 2 (n))) , 

(3.8) \\<f> - T^h ^Wh^^H^iY)) < lih (\\<j)\\L 2 (n;H 2 (Y))nL 2 (Y-H 2 (n))) 



for all (<p, 4>) S H 2 {VL) x [L 2 {Vl; H 2 {Y)) n L 2 {Y; H 2 (n)] . 

Proof. (|3.5j) and (|3.6|) are standard interpolation-error estimates, see [ID] , e.g., while 
()3.7p and (|3.8|) are interpolation-error estimates especially tailored for elliptic prob- 
lems with two-spatial scales structures; see Lemma 3.1 [B] (and its proof) for a 
statement refering to the periodic case with (n — l)-spatially separated scales. One 
of the key ideas of the proof is to see the spaces L 2 {yi,L 2 (Y)) and L 2 (fl,H 1 (Y)) 
as tensor products of the spaces L 2 (Q) and L 2 (Y), and respectively of L 2 (Q) and 
ff 1 ^). □ 

Remark 3.2. Note that, without essential differences, this study can be done in 
terms of two distinct triangulations 7~h M and 7~h m , where Km and h m are maxi- 
mum length of the sides of the corresponding triangulation of the macro and micro 
domains (51 and Y). 

Unless otherwise specified, the expressions | • | and || • || denote the L 2 and H 1 
norms, respectively, in the corresponding function spaces. 
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3.2. Main result. Proof of Theorem [3751 

Definition 3.3. (Weak solution of semi-discrete formulation) The triplet (Uq, u , v h ) 
is called weak solution of the semi-discrete formulation (|2.12[) - (|2.13p if and only if 



(3.9) / 9d t U(}(t)ipdx+ / DVU£(t)Vipdx = 

Jn 

b ((ttf + U ext - u h )(t)) d\\ + 9d t U ext {t) + DAU ext {t) ) ^dx 



(3.10) / d t u h {i)<t>dxdy + d ± V " y u N '(i)V ' y( j> dxdy = 
JnxY JnxY 

[ I b{{U% + U ext -u h )(t))<f>d\ldx-k [ r](u h (t),v h (t))(t)dydx 
JnJr R ' JnxY 

(3.11) / d t v h (t)ipdydx + d 2 V y v h (t)V v ip dydx = 
JnxY JnxY 

(u h (t),v h {t))^dydx 

JnxY 

for all ^EV h and e W h x W h and C/ O h (0) = Uj 6 L 2 {n) and u h (0),v h (0) £ 

L 2 {Q x Y). 

Lemma 3.4. (Improved regularity) Assume (Al)-(A5) to hold. Then 

(3.12) Uq £ L 2 (S ; H 2 (fl)) 

(3.13) u h ,v h £ L 2 (S; L 2 {Vt- H 2 {Y))) n L 2 (S: L 2 {Y; H 2 {n))). 

Proof. Assumption (A5) and a standard lifting regularity argument leads to Uq £ 
L 2 {S;H 2 {fl)) andu h ,v h £ L 2 {S xtt; H 2 {Y))). Employing difference quotients with 
respect to the variable x (quite similarly to the proof of Theorem 6.2 [7]), we can 
show that u h , v h € L 2 (S x Y; H 2 (fl))). We omit the proof details. □ 

Theorem 3.5. (Rate of convergence) Assume (Al)-(A5) are satisfied. If addition- 
ally, assumption (A6) holds, then it exists a constant 1C > 0, which is independent 
of h, such that 

\\Uq - U |||2( S;i? i(Q) + \\u — U h \\ 2 L2( S . L 2( n . H 2( Y )))nL 2 (S-,L 2 (Y-,H 2 (n))) 

(3-14) + \\v - v h \\\2(s-L 2 (n-H 2 (Y)))r\L 2 (s-,L 2 (Y-H 2 (n))) ^ K-h 2 . 

Remark 3.6. We will compute the constant K. explicitly; see (|3.28p . 
Proof, (of Theorem l3.5p Firstly, we denote the errors terms by 



C y 



= U -U; 



o-"o* 
u-u h 

,.h 



v — V 

We choose as test functions in Definition 13.31 the triplet 
(3.15) (<^,V) :=(r h -U h , P h -u\q h -v h ), 
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where the functions r , p , and q h will be chosen in a precise way (in terms of Riesz 
projections of the unknowns) at a later stage. We obtain 

~\U -U h \ 3 + D\\U-U h \\ 2 < [ 6d t (U -U h )(U -U h ) 
Z at Jo 



+ I D\7(Uo-US)\7(U -U h ) 

d t (U - U h )(U - r h ) + f DV(U - U h )\7(U - r h ) 
: Jn 

(3.16) - 8 I d t (U - U h ){r h - U h ) + / DV{U - U h )V(r h - U h ). 



Using Cauchy-Schwarz inequality, we have 

+ D\\U a -U h \\ 2 <9\d t (U -U h )\\U-r h \ 

+ D\V(U -U h )\\V(U -r h )\ 

+ 9\d t (U - U h )\\r h - U h \+D\V{Uv - U h )\\V{r h - U h )\ 

< 0\dt(U o - U h )\\U - r h \ + D\V(U - U h )\\V(U - r h )\ 

(3.17) + / / \b(Uo-u)-b(UH-u h )\\r h -U h \d\\. 

Jn JT R 

Noticing that r h -U = (r h - U ) + (U - U h ), (l3~T71) leads to 

d -j t \eu\ 2 + D\\ eu \\ 2 < 0\d t eu\\U - r h \ + D\Veu\\V(U - r h )\ 

(3.18) + cf f (\e v \ + \e u \) (\r h - U \ + \e v \) d\\. 

Jn Jr R 

Proceeding similarly with the remaining two equations, we get: 

^\d t e u \ 2 + dilVyeul 2 < \d t (u - u h )\\u - p h \ + dx\W(u - u h )\\V(u - p h )\ 

+ \d t (u~u h )\\p h - u h \ +d 1 \V(u~u h )\\V(p h -u h )\ 

< |$e u ||u-p h |+di|V tf e ||V w (u-p /, )| 

\b(U -u)- b(U h - u h )\\p h - u h \dXl 

in Jr R 

+ k [ \r ] (u,v)-T 1 (u h ,v h )\\p h -u h \ 
JnxY 

< \d t e u \\u-p h \ +d 1 \V y e u \\V(u-p h )\ 

+ c / (\eu\ + \e u \)(\p h -u\ + \e u \)dXl 
Jn Jr R 

+ kf \R(u)Q{v)-R(u h )Q(v h )\(\p h -u\ + \e u \). 
JnxY 

(3.19) 
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Finally, we also obtain 

/ \d t e v \ 2 + d 2 f \V y e v \ 2 < \d t e v \\v-q h \+d 2 \V y e\\\7 y (v~q h )\ 
JnxY JnxY 



(3.20) + ok \R(u)Q(v) - R(u h )Q(v h )\ (\q h - v\ + \e v \) . 

JnxY 

Putting together (|3~18l) . (|3~19)) . and ([3T20]) . we obtain 

+ di||e„|| 2 + d 2 ||e„|| 2 < 9\d t eu\\U - r h \ 

+ \d t e u \\u-p h \ + \d t e v \\v - q h \ + £»||e a |||V(C/o - r h )\ 

+ di||e„|||V0 -p h )\ + d 2 ||e„|||V,> - q h )\ 

+ c[ [ (\eu\ + \e u \)(\r h ^Uo\ + \eu\)dXl 
Jn Jtr 

(\eu\ + \e u \)(\p h -u\ + \e u \)d\l 



n Jr R 

k(l + a)\R(u)Q(v) - R(u h )Q(v h )\ (\p h -u\ + \q h - v\ + \e u \ + 



SlxY 



1=1 

where the terms Ig (£ £ {1, . . . , 4}) are given by 

h := eidteuWUo-r'^ + ldteuWu-p^ + ldte^v-q' 1 ] 

h ■= D\\Veu\\\W{U - r h )\ + d 1 \V y e u \\V y {u - p h )\ + d 2 \V y e v \\V y (v - q h )\ 

h := c [ [ (\eu\ + \e u \)(\r h -U \ + \p h -u\ + \eu\ + \e u \)d\l 
Jn Jr a 

h := k(l + a) [ \R(u)Q(v)-R(u h )Q(v h )\(\p h -u\ + \q h -v\ + \e u \ + \e v \) 
JnxY 

We choose now r h ,p h , and q h to be the respective Riesz projections of Uq , u h , and 
v h and estimate each of these 1^ terms, i.e. we set 

(3.21) r h := KffU h , p h := K%u h , and q h := K%v h . 

The main ingredients used in getting the next estimates are Young's inequality, 
the interpolation-error estimates stated in Lemma 13. 1[ the improved regularity 
estimates from Lemma 13. A\ as well as an interpolation-trace inequality (see the 
appendix in [5], e.g.). 

Let us denote for terseness 

X := L 2 (S; L 2 (n ; H 2 {Y))) n L 2 (S; L 2 (Y; H 2 (n))). 

We obtain following estimates: 
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|/i| < liO\d t e u \h 2 \\U \\ H 2 {n) +- f3 h 2 (\d t e u + d t e v )(\\u\\ x + \\v\\ x ) 

k2 7l# A~ ,2 , nrr 112 A , 1,273 /,- ,2 , 1 1 |2 



< ^ (l^ 2 + ||£/o||^ ( n)) + (IW 
(3.22) + h 2 ^{\d t e v \ 2 + \\v\\ 2 x ). 



\I 2 \ < 7 2 J5||Veu||/i||C/|| Jif 2 (n) + 7 4 di|V y e u |/»||u||x +74d2|V B e„|/i|H| x 

< e\Veu\ 2 + h 2 c el 2 D 2 \\U \\ 2 HHn) + e\V y e u \ 2 + h 2 c el 2 d 2 \\u\\ 2 x 
+ e\S7 y e v \ 2 + h 2 c € jldl\\v\\ 2 x 

< h 2 c*c e ( 7 2 + 2 7 |) (D 2 +d\+ d 2 ) {\\U \\ H , (n) + | Ml 2 , + \\v\\ 2 x ) 
(3.23) + e\Vejj\ 2 + e\V y e u \ 2 + e\\7 y e v \ 2 , 

where the constant c* > is sufficiently large. 

The estimate on |/ 3 | is a bit delicate. To get it, we repeatedly use the following 
interpolation-trace estimate 



£ 2 00)> 



(3-24) M\h m mr R ) <ef \VvV\h(Y)+c(* + l)\Mb<n: 

for the case when tp E {6 W , e v }, where e > and c,c e g]0,oo[ are fixed constants. 
We get 



|/s| < cX(T R ) I \ eu \\r h -U Q \+c I \r h -U \ \e u \d\l 



+ 



< 



I \eu\ I \p h -u\+c[ f \e u \\p h -u\d\\+2c [ f (\e u \ 2 + \e v \ 2 ) d\\ 
Jn Jv R JnJTn JnJT n 

^^(iMS^n) +71^11^011^(0)) 

+ \ (7iA(r ii )/i 4 ||J7 ||ff2 (n) + ||e u || L 2 (0;L 2 (rR)) ) 

+ \ (|A(r fl )|| ec/ || 2 H2( o) + ^MlMIx + c(c £ + i) 73 ^ 4 IMIx) 

+ 2cX(T R )\eu\ 2 + e / \\7 y e u \ 2 + c(c e + 1)|| 
(3.25) 

In order to estimate from above the term |/ 4 |, we use the structural assumption 
(A3) on the reaction terms R(-) and Q(-). We obtain 

\h\ < k(l + a) [ \R(u)-R(u h )\\Q(v)\ + \Q(v)-Q(v h )\\R(u h )\x 

X (\p h -u\ + \q h - v\ + \\e u \ + \e v \) 
< 3h 2 k l3 (l + a)(Q m c R + R m c Q ) (| Ml 2 , + Mil 2 ,) 
(3.26) + k(l + a)(Q m c R + R m c Q )(\e u \ 2 + \e v \ 2 ) , 
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where R m '■= max, re [ M2 ]{R(r)} 7 Q m := max se [ 0! M 3 ]{-R(s)} > while cr and cq are 
the corresponding Lipschitz constants of R and Q. 
Consequently, we obtain 



< h 



dm 



'teu 



73 



*±\d t e u Y + \d t e v Y)+h*{1C + T{h)) 



k(l + a)(Q m c R + R m c Q ) + c + 



(3.27) 
where 



+ e|Vec/| 2 



e 2 



|V«e v | 



/C := ^ 



f (Nli 



(3.28) 

h 2 

(3.29) 



Ti 



\v\\ 2 x ) 



-^-H t/ o||^2 (n) + 3fc7 3 (l + a{Q m c R + R m c Q ) (\\u\\ 2 x 



\v\?x) 



2e74||u||i+ 7 i 1 + -A(r fl ) \\U { 



y ollff 2 (n) 



2c(c e + l)7 3 ||M| 



A" • 



Notice that K. is a finite positive constant that is independent of /i, while F :]0, oo[— > 
]0, oo[ is a function of order of 0(h 2 ) as h — > 0. 

By (A6) we can compensate the first term of the r.h.s. of (|3.27[) . while the 
last three terms from the r.h.s. can be compensated by choosing the value of e as 
0, min{Z), d 2 > "^4} ■ Relying on the way we approximate the initial data, we 



e G 



use now Gronwall's inequality to conclude the proof of this theorem. 



□ 
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